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We consider general relativistic Cauchy data representing two nonspinning, equal-mass black holes boosted toward each other. 
When the black holes are close enough to each other and their momentum is sufficiently high, an encompassing apparent horizon 
is present so the system can be viewed as a single, perturbed black hole. We employ gauge-invariant perturbation theory, and 
integrate the Zerilli equation to analyze these time-asymmetric data sets and compute gravitational wave forms and emitted 
energies. When coupled with a simple Newtonian analysis of the infall trajectory, we find striking agreement between the 
perturbation calculation of emitted energies and the results of fully general relativistic numerical simulations of time-symmetric 
initial data. 
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I. INTRODUCTION 

The collision of two black holes is expected to be an im- 
portant source of gravitational radiation for gravity wave 
detectors currently under construction. A major theoret- 
ical effort is underway to compute the gravitational wave 
form from the orbiting inspiral and coalescence of black- 
hole binaries. Post-Newtonian theory should provide a 
sufficiently accurate wave form for much of the inspiral 
phase and enable the extraction of considerable informa- 
tion about the parameters of the system jL|. However, it 
is anticipated that a fully general relativistic treatment 
will be necessary to predict the wave form from the final 
stages of coalescence. 

Recently, the seminal calculations of Smarr and Epp- 
ley Q for the head-on collision of two equal-mass holes 
starting at rest from a finite separation have been redone 
H with the benefit of new numerical techniques and the- 
oretical tools, as well as vastly increased computational 
resources. The basic conclusions reached by the modern 
calculations are remarkable for their similarity to those 
of the original study. The maximum amount of energy 
radiated is small, less than 0.1% of the mass of the sys- 
tem. Also, the wave forms are indistinguishable at the 
level of the numerical accuracy from quasinormal mode 
oscillations of a black hole. 

Compelled by these results, Price and Pullin Q an- 
alyzed the Misner initial data for two black holes at a 
moment of time-symmetry as if it represented a single, 
perturbed black hole. Using gauge-invariant perturba- 
tion theory and the Zerilli equation, they were able to 
compute the initial distortion of the black hole, and the 
resulting asymptotic wave forms and energy fluxes as a 
function of separation. For small separation, when the 
approximation of the merged system as a single, per- 
turbed black hole is expected to be most valid, the agree- 
ment between the radiated energy from perturbation the- 
ory and the results of fully relativistic simulations Q is 



excellent, apparently within the error bars of the numer- 
ical calculation. The quadrupole wave forms are also re- 
markably similar when read off at the same radius. Only 
when the initial separation of the black holes is somewhat 
larger than the cutoff for encompassing apparent horizons 
is there substantial discrepancy between the perturbation 
theory and evolution results for radiation efficiency. 

Previously (Ref. ||, hereafter Paper 1), the current au- 
thors had used a similar perturbation theory analysis to 
study spurious radiation in time-symmetric and asym- 
metric two-black-hole data sets in the opposite limit - 
that of large separation. These data sets represent the 
direct extension of Misner data allowing the two black 
holes to have nonvanishing initial linear and angular mo- 
menta. In this paper, we apply these gauge-invariant per- 
turbation techniques, in the spirit of the Price and Pullin 
paper, to the close limit when the two black holes have an 
encompassing apparent horizon, and examine the gravi- 
tational radiation wave form and the amount of energy 
radiated for time-asymmetric initial data. 

There are several motivations for this study of boosted 
black-hole initial data. In the final plunge phase of bi- 
nary black hole coalescence, the black holes are likely 
to have substantial infall velocities. Although the ac- 
tual collision is not expected to be head-on, axisymmetric 
calculations of boosted head-on collisions are an interest- 
ing limiting case of full three-dimensional plunge simula- 
tions. In addition, these calculations extend the regime 
where perturbation calculations can be fairly compared 
with the fully relativistic simulations. It is clear that, for 
time-symmetric initial data with no event horizon, the 
perturbation theory calculation should greatly overesti- 
mate the distortion of the merged black hole, and thus 
the radiated energy. For these cases the merged black 
hole doesn't form until the individual black holes have 
evolved toward each other and have some infall momen- 
tum. It is this merged black hole that one would like to 
analyze with perturbation theory and compare against 
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the evolution results. 



II. METHODOLOGY 

In this paper, we will be concerned with two equal- 
mass nonrotating black holes, each with axisymmetric 
inward-pointing momentum P (the slice has zero net- 
momentum). Initial-data sets representing one or more 
black holes with individually specifiable linear and an- 
gular momenta are constructed using the conformal- 
imaging approach developed by York and coworkers 
The case of two black holes with axisymmetric momenta 
was implemented numerically by Cook JtJ . The data sets 
used for this study were constructed using a code based 
on the Cadez-coordinate approach described in that 
work. The reader should refer there for details regarding 
the construction of the data sets and for further descrip- 
tions of the parameterization of the data sets described 
below. The separation of the holes is parameterized by 
(3, which is related to the bispherical-coordinate separa- 
tion parameter /j,q by the relation hq = cosh _1 (/3/2) in 
the case of equal-mass holes. The code is used to com- 
pute the inversion-symmetric (with minus isometry con- 
dition) extrinsic curvature Kij and to solve the Hamil- 
tonian constraint for the conformal factor ip. Once the 
full initial data is computed, several physical quantities 
characterizing the system are computed. Of interest in 
this paper are the ADM mass of the initial-data slice M, 
the proper separation of the holes £, and the masses of 
the individual holes mi = TO2 defined in terms of the 
area of the marginally outer-trapped surface associated 
with each hole. We also define the total or bare mass of 
the system by m — mi + m-i- Note that the difference 
between m and M is due to the binding energy of the 
system. Given an initial-data set, we use the boundary- 
value-problem method described in Paper 1 to locate all 
marginally outer-trapped surfaces surrounding the two 
holes (if they exist) and identify the apparent horizon(s). 
We should note that there is nothing unique about our 
choice of initial data for representing two colliding black 
holes. One could imagine, for example, considering data 
with Euclidean topology with the black holes represented 
by boosted matter collapsed inside its horizon. Our ini- 
tial data was chosen for the convenience of its highly 
refined numerical treatment and earlier physical ex- 
ploration (Paper 1). 

Like Price and Pullin Q , for purposes of our analysis 
we treat the spacetime as a perturbation (but not neces- 
sarily a time-symmetric one) of Schwarzschild. First, we 
establish a Schwarzschild-like coordinate system around 
the two black holes in terms of the (background space) 
isotropic coordinates used in the numerical solution r = 
ri(l + M/2ri) 2 . Note that the background space of the 
numerical solution can be directly parameterized by the 
isotropic radial coordinate ri even thought the numer- 
ical solution is found in Cadez coordinates. The total 



ADM mass of the slice M is used as the Schwarzschild 
background mass. Tortoise coordinates are also defined 
in the usual way: r* = r + 2M \n(r /2M — 1). Com- 
putation of wave perturbations involves the calculation 
of multipole amplitudes by surface integrals. These are 
performed over constant Schwarzschild radial-coordinate 
two-spheres. The integrands involve the conformal fac- 
tor ip, Schwarzschild-coordinate extrinsic-curvature com- 
ponents K^, and their Schwarzschild-coordinate radial 
derivatives. Calculation of these quantities at their re- 
quired locations is achieved with bi-cubic spline interpo- 
lations and a series of coordinate transformations. 

The gauge-invariant function Qi m is formed out of 
multipole projections of ip and ip^ computed by nu- 
merical integrations over a coordinate two-sphere (cf. 
Refs. |PJs|JTo[|). For this paper we compute only the case 
of Q = Q20- We also require the Schwarzschild time- 
derivative of the gauge- invariant function dtQ. This time 
derivative is computed as 

d t Q = aC n Q, (1) 

where £ n is the Lie derivative along the slice-normal con- 
gruence n and the factor a — — 2M/r) comes about 
from the transformation from the slice-normal time co- 
ordinate to the Schwarzschild time-coordinate. The Lie 
derivative of Q is calculated using the extrinsic curvature 
(and its radial derivative) via the definition 

C n9ll = -2K ir (2) 

The gauge-invariant perturbation function and its time 
derivative, known as a function of radius surrounding the 
merged black hole, serves as initial data for integration 
of the Zerilli equation. The numerically generated initial 
perturbation is interpolated onto a fine grid (typically 
8000-16000 zones) that is even in r* and extends from 
n = -500Af to r* = 2000M. The Zerilli equation (cf. 
Ref . |]Jl(],^| ) is then integrated forward in time until the 
whole perturbation has been propagated to |r*| — > oo. 
Approximate asymptotic wave forms and energy fluxes 
are computed at large radii. 

Our code for calculating the initial black-hole per- 
turbation from numerically generated initial data was 
checked by comparing it against the time-symmetric re- 
sults of Price and Pullin Q . It should be noted that they 
analytically expanded the metric perturbation about 
Schwarzschild in powers of the parameter e = 1/] In yuo | 
and retain only the leading term in e. In the limit of 
small separation, the initial perturbation we obtain nu- 
merically agrees closely with their analytic results (for 
(3 < 3.25 the agreement is better than 5%). Not surpris- 
ingly, for larger separations the differences become larger. 
For the horizon cutoff point of (3 = 4.17 or /iq = 1.36 
(the largest separation allowing an encompassing appar- 
ent horizon for time-symmetric initial data), the analytic 
prediction for the energy radiated is about 60% higher 
than the result from the full solution. It is interesting to 
note that neglecting the higher-order terms in e always 
seems to lead to a greater amount of radiated energy. 
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FIG. 1. The apparent horizon formation line. The inward 
linear momentum on each hole P/m is plotted as a function 
of proper separation i/m. 

III. WAVEFORMS AND ENERGY FLUX 

The perturbed black-hole approximation assumed in 
this paper is not valid if the two black holes have not 
merged (have no common event horizon). For separa- 
tions small enough that an encompassing apparent hori- 
zon exists, we find that the addition of inward-pointing 
linear momenta makes the encompassing apparent hori- 
zon more spherical and the maximum radiation efficiency 
(defined as the ratio of M minus the mass of the apparent 
horizon to M) decreases. Moreover, we find that the met- 
ric perturbation Q always gets smaller. We, therefore, 
contend that our treatment of time-asymmetric initial- 
data sets with inward-pointing linear momenta is always 
at least as valid as the study of the time-symmetric so- 
lution. 

In Paper 1 we located the horizon-formation line for 
boosted black-hole initial data. For a given separation 
parameter (3, we searched for the smallest value of in- 
ward linear-momentum for which an encompassing hori- 
zon surrounded the holes. The momentum as a function 
of proper hole separation for this horizon line is displayed 
in Fig. [l]. The horizon cutoff point mentioned previously, 
at — 4.17 and P = 0, lies on this line. For larger values 
of /3, it is necessary to give the holes inward momentum in 
order for an encompassing apparent horizon to exist. We 
note that this horizon line is only an estimate of where the 
actual encompassing event horizons will form. Along this 
horizon-formation line, we have computed the radiated 
energy for the initial-data sets using the gauge-invariant 
perturbation formalism and Zerilli-equation integration 
method described above. In Fig. ||, the radiated energy 



FIG. 2. Radiated energy along the apparent horizon forma- 
tion line. The radiation efficiency (the total radiated energy 
as a fraction of the ADM mass) computed with perturbation 
theory is plotted for initial data along the horizon line, pa- 
rameterized by the proper separation of the holes £/m. In 
the inset we show the radiation efficiency plotted as a func- 
tion of the separation of the corresponding time-symmetric 
initial-data set, (£/m)o, computed using Eq. (3). 

is plotted as a function of proper separation £/m. For the 
points shown, the inward momentum on each hole ranges 
from P/m — 0.0355 to P/m = 1.738. One interesting fea- 
ture is that the radiation efficiency appears to saturate 
at about 2%, substantially below the maximum radia- 
tion efficiency based on area theorem arguments. This 
suggests that it may be impossible to obtain high radia- 
tion efficiency for black-hole collisions, even if they merge 
with very large momenta. As a gauge of what constitutes 
a large momentum, we estimate below the momenta of 
two holes at the moment of horizon formation assuming 
a parabolic infall from rest at infinite separation. 

At each point on the horizon formation line, we have 
a value for the separation t/m and momenta P/m of the 
holes. Treating the black holes as point particles and us- 
ing Newtonian dynamics, we can estimate the separation 
(£/m)o at which the holes were at rest: 



\mJ Q 1 - 8{P/m) 2 (£/m)' y ' 

Clearly, if we assume infall from rest, the maximum mo- 
menta the two holes can obtain at the point of hori- 
zon formation is estimated by the point on the horizon- 
formation line where the denominator of Eq. (0) vanishes, 
i.e. (£/m) — > oo. We find this point to be P/m — 0.249 
and i/m = 2.01. At this point, the implied radiation effi- 
ciency is less than 0.16%. In the inset of Fig. ||, we show 
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form. The presence of momentum (extrinsic curvature) 
on the initial slice causes a significant transient feature 
in the wave form; this qualitative effect is not seen by the 
numerical relativity simulations of time-symmetric initial 
data with fj,o = 2.2 because it takes finite time for the mo- 
mentum component of the perturbation to propagate out 
to the extraction surface. It is possible that this feature 
will be present in the evolution of time-symmetric initial 
data starting at larger separations. 



FIG. 3. Waveforms from time-symmetric and asymmet- 
ric two-black-hole initial data. The top curve shows the 
quadrupole wave form from analysis of a time-symmetric data 
set with (3 = 4.275. The lower curve shows the quadrupole 
wave form from analysis of a time-asymmetric initial data set 
with the same value of j3 and P/m = 0.1225. Both wave forms 
are extracted at a radius r = 200M. 

the radiated energy from points on the horizon-formation 
line plotted as a function of the Newtonian estimate for 
their proper separation when at rest. We find striking 
agreement between our calculation of total energy radi- 
ated and the simulations of Anninos et al. || . For cases 
where their initial data had an encompassing horizon, it 
is clearly correct to compare with the perturbation anal- 
ysis of Misner data. For cases with greater separation, 
the time-asymmetric analysis gives excellent results. For 
example, for [1q = 2.2 or (£/m)o = 3.97, the radiation 
efficiency from perturbation theory of the corresponding 
horizon- line initial-data set is E/M = 7.9 x 10 -4 , as com- 
pared with E/M = 1.7 x 10~ 3 from the time-symmetric 
analysis of = 2.2 Misner data and E/M ~ 5.5 x 1CT 4 
from the fully relativistic simulations. A post-Newtonian 
calculation of the infall trajectory might improve this 
comparison. 

In Fig. ^ we show a typical wave form from a boosted 
head-on collision observed at a radius of r = 200 M. The 
case shown is for a separation of /3 = 4.275, P/m = 
0.1225, and i/m = 1.948. For comparison we show the 
wave form from time-symmetric data with the same (3. 
Both wave forms are dominated by normal mode oscil- 
lations within about 10M after the black-hole surface is 
causally apparent at the extraction radius. The addition 
of ingoing momentum considerably increases the ampli- 
tude of the oscillation and reverses the sign of the wave 



IV. DISCUSSION 

Anninos et al. |J have shown that for collisions result- 
ing from large initial separations there is excellent agree- 
ment between the emitted energy and the well-known 
results for a test-particle falling into Schwarzschild cor- 
rected for equal-mass objects, finite infall distance, and 
horizon heating. Price and Pullin Q demonstrated that 
a perturbation analysis of time-symmetric initial data 
could reproduce the results of the fully relativistic simu- 
lation in the case that the black holes have small initial 
separation. Here we have shown that the perturbation 
analysis can be extended to larger separations, includ- 
ing the regime in which point particle analysis is valid, 
by considering appropriate time-asymmetric initial-data 
sets. Adopting this perspective, one can accurately pre- 
dict the total emitted energies over the entire range of 
separation, from the close limit to parabolic trajectories 
starting at infinity. 
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